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Abstract. We illustrate an efficient new method for handling polynomial systems 
with degenerate solution sets. In particular, a corollary of our techniques is a new 
^ • algorithm to find an isolated point in every excess component of the zero set (over 

pLn , an algebraically closed field) of any n by n system of polynomial equations. Since 

QQ ' we use the sparse resultant, we thus obtain complexity bounds (for converting any 

input polynomial system into a multilinear factorization problem) which are close 
to cubic in the degree of the underlying variety — significantly better than previous 
bounds which were pseudo-polynomial in the classical Bezout bound. By carefully 



o 



•^ ■ taking into account the underlying toric geometry, we are also able to improve 



^ , the reliability of certain sparse resultant based algorithms for polynomial system 

solving. 



1. Introduction 



(N 



The rebirth of resultants, especially through the toricQ resultant 
^ ■ begun to provide a much needed alternative to Grobner basis methods for solving 

C^ ! polynomial systems. Continuing this philosophy, we will present a new, fast and 

reliable, resultant-based method for handling certain degenerate polynomial systems. 
GN ' Simply put, we refine and generalize the generalized characteristic polynomial 

(GCP) [|Can9CI| , |Shu93|] to take sparsity into account. Furthermore, we introduce 



j^ ! the concept of a twisted Chow form in order to completely avoid any degeneracies 

S I within our algorithm. 

The importance of dealing with degenerate polynomial systems has been observed 
in earlier work on quantifier elimination over algebraically closed fields ||CG84| , [Ren87 



r^ I |Can88|| : Many reasonable algorithms fail catastrophically when presented with an n 

d [ by n system having positive-dimensional zero set. When such is the case, it is of 

considerable benefit to the user to at least be given some sort of description of the 
zero-dimensional part of the zero set. This was a benefit of Canny's original GCP, 
but he remarked [|Uan9CI| , pg. 242] "For large, dense problems however, the resultant 



and GCP methods should be faster [than Grobner algorithms]." Our construction, 
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the toric GCP, promises to be much more competitive when apphed to sparse 
systems in such a comparison. 

Remark 1. It should be emphasized that perturbation methods for degenerate sys- 
tems (such as the toric GCP) are of the greatest importance when working with ex- 
act arithmetic. However, floating point polynomial system solving also benefits from 
a complete and rigourous understanding of the potential degeneracies within exact 
arithmetic, e.g., |^ma87|| . In any case, it is frequently the case in many applications 
that "real life" happens to land in a measure exception which breaks an algorithm. 

In what follows, we will frequently use multi-index notation in order to precisely 
control which monomial terms are allowed to appear in our polynomial systems.0 This 
notation, and in particular, supports and Newton polytopes, are amply detailed 
in earlier works of Emiris |[Emi94|| , Huber ||Hub96|| , Gelfand et. al. ||GKZ90| , |GKZ94 



Rojas ||Roj94| , |Roj97b|| , Sturmfels ||Stu93|| , and Verschelde || Ver96|| , to mention but a 
few authors and references. 

Since what we will present is at heart a perturbation method, we will first need the 
following definition to construct certain polynomial systems in "general position." 



Definition 1. |[Roj94| , |KW^ ] Given n-tuples D -.= {01, .. . , D„) and E:={Ei, . . . , En 
of nonempty finite subsets of Z", we say that D fills E iff (0) Di C Ei for all i G [n] 
and (1) M.{D) = M.{E). An irreducible fill is then simply a fill which is minimal 
with respect to n-tuple containment. If D or E is instead an n-tuple of polytopes in 
W^ , then we will use the same definition. 



In the above, M{E) denotes the mixed volume pZSBj §cEg, (EU^^I , PGH96|| of 
the convex hulls of the Ei, and we use [j] for the set of integers {!,... , j}. 

Our construction is summarized in the following definition and main result. (Hence- 
forth, all of our polynomials and roots are to be considered over an algebraically closed 
field K.) 

Definition 2. Suppose F is an nxn polynomial system with support contained in 
E and g{x) '■=J2e€A'^e^'^' ^here the u^ are algebraically independent indeterminates 
and AgIj"' is nonempty and finite. Assume further that M.{E) >0 and A has at least 
two elements. Letting D be an irreducible fill of E, F* := (Xles-D -'"'^ I ^ ^ W)' '^'^'^ 

Res^E,A){F, g) and H{u; s) ■.= Res^E,A)iF - sF*, g). 
We call Ti a toric generalized characteristic 



u:={ue I e&A), define C\ia{u):= 
where s is a new indeterminate. 
polynomial for (F,A). 



In the above, Res^(-) denotes the toric resultant. Recall also that to any ?7,-dimensional 
rational polytope P C M" one can associate its corresponding toric variety Tp 
PiE7g| , PEZ^ , F;IK)3| , [Roj97b|| . 



^This is what we mean by sparsity, as opposed to working solely in terms of polynomial degrees. 
Note that our notion overlaps the case of having very few monomial terms. 
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Main Theorem. Following the notation of definition ^ consider Ti as a polynomial 
in s with coefficients in K[u\ and let ^a{u) be the coefficient of the lowest degree 
term ofH. Also let F:= ^^^^ Conv(^i), P := P + Conv{A), and ^ : Tp — > pl^'"^ 
the natural morphism defined by x \-^ [x^ \ e G A] (cf. section ^. Then Ta is a 
homogeneous polynomial, of degree Ai{E), with the following properties: 

0. The constant term of7i{s) is precisely Ch^. In particular, letting Z be the zero 
set of F in Tp, Chyi = <^==^ '^{Z) is positive- dimensional. 

1. If C & Tp is an isolated root of F then Ta is divisible by XleeA'^e'^e? where 
[ce\eeA] = ^{C). 

2. If a nonzero linear form XlegA'^e'^e divides Ta then [ce | eGv4] = y9((^) /or some 
rootCeTp of F. 

3. Ta splits completely into linear factors. In particular, under the correspon- 
dence of (2), we can explicitly find at least one root of F within every positive- 
dimensional irreducible component of the zero set of F inTp. 

We call T A a toric perturbation of Ch. a- 

A complete criterion for finding fills is given in theorem |I| of section ^ and some 
simple examples of filling (and our Main Theorem) appear in section |^. Compatibility 
is defined in section ^. There we also detail what is meant by the zero set of F in a 
toric variety. 

Remark 2. Alternatively, if one wants to avoid filling, one can substitute forF* any 
polynomial system (with support contained in E) known to have exactly M.{E) roots 
(counting multiplicities) in (K*)"' . So one can also view filling as a deterministic way 
of explicitly constructing a "generic" polynomial system. 

So how does one actually use the above theorem? One simple example is the 
sparse u-resultant | |Emi94] ] , which is simply a variant of the classical w-resultant 



|Van50| | . It can be defined simply as Ch^(w) where A is the vertex set of the standard 
n-simplex in M". One useful (and easily verified) property of the sparse u-resultant is 
that F has a root (:= (Ci, ... , Cn) ^ (iT*)" =^ Uq + CiWi + ■ — h CnUn divides Ch^(w). 
So, assuming ChA{u) is not identically (and that one has good software for toric 
resultants and multivariate factoring), one can find the isolated roots of F simply by 
factoring Ch a{u). Thus the degenerate instances Chyi(w)=0 obstruct this reduction 
to factoring and our Main Theorem allows us to avoid this problem: Simply use 
Ta{u) instead of Ch a{u). 

Remark 3. An interesting 'failure" for the classical u-resultant is the case where F 
has only finitely many roots in affine space, while having infinitely many in projective 
space. This is where the classical GCP is especially handy. For the toric GCP, the 
toric variety Tp plays the role of projective space. This is part of the philosophy of 
sparse elimination theory: working in a well chosen toric compactification (depending 
on F) leads to better algorithms than if one were to work only in projective space. 
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In particular, lifting to Tp helps us detect precisely when the sparse u-resultant is 
identicaly 0. 

However, there is a stickier subtlety which occurs with the sparse -u-resultant: It is 
possible for F to have only finitely many roots within Tp with Ch^(ii) still vanishing 
identically (cf. section |2.3|) . This is the motivation for twisted Chow forms, which 
are defined below, and pursued at much greater length in ||Roj97c |. 

Corollary and Definition 3 (cf section ^. Following the notation of our Main 
Theorem, pick A to he the vertices of a product of simplices with which P is com- 
patible. We then call Chyi(n) the twisted Chow form of the zero set of F with 
respect to A. Furthermore, a twisted Chow form does not vanish identically if F 
has only finitely many roots in Tp. D 

In particular, since our last construction implies that '7conv{A) is a product of 
twisted projective spaces |[Ful93|| , the coefficients of a twisted Chow form are ac- 



tually multisymmetric functions [pS95| , [Roi97a|l of projections (by (/?) of roots of 
F (in Tp) onto this product. So twisted Chow forms generalize the w-resultant, 
the sparse ti-resultant, and (suitably extended |[Roj97c|| ) the Chow form of a pro- 



jective variety |PS95|| . Moreover, when combined with the Smith normal form 
[111891, [HM91| , [HS95II and the toric GCP, twisted Chow forms give a considerably more 
reliable method for solving polynomial systems than the sparse -u-resultant |[Roj97c|| . 
An important difference to note is that our present toric GCP is primarily suited 
for finding roots in [K*)^, while the original GCP is mainly suited (in a "non-sparse" 
way) for affine space. To completely generalize and improve the GCP in affine 
space, it is necessary to use the aflflne sparse resultant and this is pursued further 



m 



For instance, by replacing the sparse resultant with the affine sparse 
resultant, and using i^^-counting [|lloj97b|| instead of ffUing, we can actually recover 
Canny's GCP in the dense case. 

We close this introduction with a word on the computational complexity of com- 
puting the toric GCP. Neglecting preprocessing (ffnding an irreducible fill and finding 



a mixed subdivision in order to set up the toric resultant matrices | |EC95|| ), recent 



work of Emiris, Morrain, and Pan |[EP97], |MP97|| suggests that it is possible to find 



the perturbation J-'a within a number of arithmetic stepsQ which is close to cubic 
in Ai{E). (Indeed, Canny has pointed out | ]Can9U[ | that the original GCP can be 
computed in time close to cubic in the Bezout bound, which would be a special case 
of time cubic in A4{E).) Since Ai{E) is much smaller than the Bezout number for 
most polynomial systems [[Roj94| , [Roj97b|| , these preliminary results suggest that the 



toric GCP has considerable potential for practical applications. 
Let us now illustrate some of our theory. 

^Althou gh the possibi lity of "near-quadratic" algorithms for polynomial system solving is dis- 
cussed in [EP97, |MP97 |, neither paper discusses the case where the zero set of F is positive- 



dimensional, which is our main concern here. 
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2. Examples 

We begin with two small examples of filling. We then see applications of the toric 
GCP to some degenerate 2x2 and 3x3 polynomial systems. Finally, we see a brief 
comparison of the toric GCP to the original GCP. 

2.1. Filling Squares and Cubes. For our first example, consider the pair of rect- 
angles P:= ([0, a] X [0, 6], [0, c] X [0, rf]) where a, b, c, and d are positive integers. Then 
it is easily verified (via theorem |I] of section ^ that the pair 

D:=({(0,0),(a,6)},{(0,d),(c,0)}) 

fills V. In this case, the mixed area of both pairs is easily checked to be ad + bc. Note 
also that D is a pair of oppositely slanting diagonals of our initial pair of rectangles 
(modulo taking convex hulls). Finally, it is easily checked that D is indeed irreducible, 
since the removal of any point of D results in a mixed area of 0. 

For our second example, let V instead be a triple of standard cubes (so that the 
vertex set of each cube is simply {0, 1}^). Then, using the criterion from theorem |l| 
once again, it is easily verified that the triple 

D:=({(0, 0,0), (1,1,1)}, {(1,0,0), (0,1,0), (0,0,1)}, {(1,1,0), (1,0,1), (0,1,1)}) 

fills V. Also, it is easily checked that the mixed volume of both triples is 6. Note that 
the triple D consists of a body diagonal and two oppositely oriented (but parallel) 
sub-triangles of the unit cube, modulo taking convex hulls. Finally, note that this D 
is irreducible as well by theorem |I|. (This is also easily checked by one of the publically 
available software packages for mixed volume by Emiris, Huber, or Verschelde.) 

2.2. A Degenerate 2x2 System. Consider the bivariate polynomial system F: = 
(1 -|- 2x — 2x'^y — 5xy + x"^ + Sx^y, 2 + Qx — 6x^y — llxy + Ax^ + 5x^y). Letting E 
be the support of F, the reader can easily verify that Ai{E) = 4, and that the only 
roots of F are the points {(1, 1), (|, |)} and the line {—l}xK (assuming charA'7^2).0 
So it would appear that the w-resultant (and even the sparse w-resultant) will vanish 
identically and not give us any useful information about any of these roots. Let us 
see if perturbing the sparse n-resultant helps... 

Note that by theorem [|, D= ({O, (3, 1)}, {(1, 1), (2, 0)}) is an irreducible fill of E. 
So applying our Main Theorem with this D (and A:={{0,0), (1,0), (0,1)}) one can 
compute with the use of Maple that H{u] s) is precisely 

(U2 — U(j + Ui + bUiU2 — 4ltlM2 — 4UiM2JS 

+ (— 20m2'"() — 20u2U() — 4uiuq + 36uiU2 — 19u(, — 24m2 + (!iUqUiU2 

o A 2 2222 2 3 3\T 

+36ttiit2 + 36ui — 12uoUiU2 — QmiUq + 3m2'"() + 'i&uou\U2 — Auou-i — 84uiU2)s 
+ (-170u2Mo - 394u?U2 - 98uiUo - 98uo"iU2 - 20uo + 220^2 + ZlQuluo 



*When charier = 2, the second isolated root becomes an isolated root lying on the x-axis. 
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+ 14uouiU2 — llOuoUi — 226uiUq — 354miM2 + 454ui — 27AuquiU2 + 74miM2)s 

+ {-l&l2uluiU2 + 1008-U2MO + 903?io - 624uiUo - 2&'i2uluo - 2lMuoulu2 - 970^2 

-lOlOuiui + 418u?U2 - 2104uoUiMi - 642u?ui - 1547u?uo - 936uou? - Ihhlulul + 2204mi)s^ 

+{538uoMiU2 + 1271^0 + 12253^2 "o + 6972u2?io + 1929wi - 3075^1^2 + QbAuouml 

+50uiMo + 2156U2 - 960^1^0 - 2290mom? + 132uiU2 - 5344uom?W2 - lU2u\u2 + 8708^2 wo)s* 

+(4384uiuo - 24988uiMo - 1582u?u2 - 6756^0 + IQ^Muomul + 3802^1^2 + 15438uott?w2 

+ 1024«oM? + 8324u?iio - 12826Miuo + 11270uo«iU2 - 6976ut + llUulul - 21326u2ug - 2408u2)s^ 

+(3436uf U2 + 3800U0W? + 7756^1^0 - 3886^1^1 + 1225u2 + 17059uiuo - 5984u?uo 

+15708u2Mo - 12232«oWiit2 + SISOuq - 2miu\ul - 6828mom?W2 + l?,l6u\ - 12700itoWiU2 - 4312itiUo)s^ 

+ (-196^^ - U^ul - 768ulu2 + bUuful - 1288uluo 
—2A36u2Uq + I92OU0U1U2 + 1536uqUiU2 + 102AuoUiU2 

+384uiul - 1792u2ul - 38Auoul + 768uiul + 260ulul - 64iz^)s. 

(Sparse resultants of this size are quite amenable with the aid of Maple, following 
the technique of a similar computation in ||Koj97a| .) So our JF4 is just the coefficient 



of s or s^ in this polynomial, according as charA'7^2 or charii' = 2. To simplify our 
discussion, let us henceforth assume the former possibility. 

Factoring with Maple, we obtain that JF4 can be written as follows: 

-4(mo + Ml + -U2)(28mo + 4^1 + 49^2) (wo - "Ui + M2)(4mo - 4^1 + ^2) 

In particular, given any factor above, the ratio of the coefficients of Ui and Uq is 
precisely the i— coordinate of a corresponding root of F. Thus the first two factors 
correspond precisely to the two isolated roots we already know. As for the last two 
factors, note that they both give isolated points lying on the aforementioned line 
{—l}xK. This can be interpreted as assigning an excess intersection multiplic- 
ity of 2 to the line, so that the sum of all intersection numbers (of the irreducible 
components of the zero set of F in Tp) is 4. 

Note that the original GCP could have been used above but would have resulted in a 
ti-form of degree 16 (the product of the degrees of /i and /2). Also, the corresponding 
version of 7i(-) is significantly larger, having 672 terms, compared to 110 for our above 
toric GCP. 
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2.3. Stranger Degeneracies. Here we give two examples showing how the sparse 
u-resultant can fail to find roots in Tp, even with the benefit of the toric GCP, unless 
some other construction (such as a twisted Chow form) is used. 

First consider the parameterized bivariate system F := {aiy + a2X + a^xy, biy + 
b2X + b^xy). Note that the mixed volume bound for this system is 1. The sparse 
u-resultant for this system is also easily found (using the same techniques as in our 
last example) to be: 

(0362^1 + bla2ai - 63020361 - bj,a^b2ai)uQ 
+ (6^0203 - 6i63a2a.i - 62010361 + b2b^a\)ui 

+ (616302 - 62020361 + 620301 - 6263020i)m2 

In particular, we see that when 

(oi, 02, 03, 61, 62, 63) = (0, 1, 2, 0, 1, 3). 

the sparse w-resultant vanishes identically. For this specialization, it is also easy to see 
that F has only one root in 7p, and this root lies at the point of Tp corresponding to 
the vertex (0, 1) of P (cf. lemma |T] from section ^, following the notation of our Main 
Theorem. In fact, for this P, Tp = ¥\ and the vertex (0, 1) corresponds precisely to 
the point [x : y : 2] = [1 : : 0] in this particular copy of ¥\ {z denoting an extra 
variable for homogenizing). 

There are two ways of viewing this degeneracy {Ch.A{u) = while F has only finitely 
many roots in Tp) of the sparse w-resultant. The first is pragmatic: One should not 
count this as a deficiency of the sparse w-resultant because we've cheated and set two 
of the coefficients of F to 0, thus changing the supports. (Our next example avoids 
this trick.) The second point of view is more geometric: By our Main Theorem, this 
particular sparse resultant must vanish due to the fact that our specialization of F 
results in the existence of infinitely many roots of F in 7^. In fact, F vanishes on 
the 1-dimensional subvariety of Tp corresponding to the left-hand vertical edge of the 
hexagon P (cf. section^). 

Going one dimension higher, consider instead the 3x3 system G, consisting of the 
following polynomials: 



01^2; + a2xz + a^xy + a^^xyz 
biyz + 62x2; + 632;^ + b^xyz 
ciyz + C2XZ + c^xy + c^xyz 



Note that the mixed volume bound for this system is again 1. 

Clearly, -^G is a linear system in {-, -^ -}. So by Cramer's rule, we can express x, 
y, and z as ratios of 3 x 3 determinants in the coefficients. Combining this with the 
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product formula for toric resultants |PS93|| (and clearing denominators) we obtain 
that the sparse w-resultant of G is precisely^ 

[423][143][124]mo + [123][143][124]mi + [123][423][124]m2 + [123][423][143]m3 
where the bracket [ijk] [ pS95|| is the 3x3 sub determinant 



tti 


Uj 


ak 


k 


b, 


bk 


Ci 


C] 


Ck 



of the coefficient matrix of G. This compactly expressed resultant can be thought 
of as a semi-mixed Chow form — a toric resultant of a semi-mixed system [[HS95 
compressed in terms of suitable brackets. 
Now consider the specialization of G to 

yz + xz + 2xy + 3xyz 

yz + xz + Axy + S^xyz 

yz + xz + 8xy + 27xyz 



It is then easily verified that G has exactly one root in 7p = P|.: 

[x:y:z:w] = [l:-l:0:0] 

{w denoting an extra variable for homogenizing) .0 More to the point, the sparse 
w-resultant vanishes identically for this specialization of G, even though G has no 
zero coefficients. Furthermore, one can easily check that the correspondence of (2) 
(from our Main Theorem) does not give us the root of G in Tp. (Using D : = 
({(0, 1, 1), (1, 1, 1)}, {(1, 0, 1), (1, 1, 1)}, {(1, 1,0), (1, 1, 1)}), the coefficients of the u^ 
in J^A suggest the nonsensical root [0 : : 5 : 21].) 

To remedy this, we can use the twisted Chow form Ch.A'{u) with 

A':={(0,1,1), (1,0,1), (1,1,0), (1,1,1)}. 

In particular, when the coefficients of G are unspecialized. 



ChA'(u) 



Cl\ Ci2 Cij, CI4 

bi 62 63 64 

Ci C2 C3 C4 

""(0,1,1) ^(1,0,1) ""(1,1,0) ^(1,1,1) 



^We also need the fact that the Pedersen-Sturmfcls formula, o riginally stated only over C, remains 
true over a general algebraically closed field. This is proved in | Roj97d |. 



^ If charii" S {2,3} then G will actually have infinitely many roots in Tp. So let us assume 
henceforth that chari^ ^ {2, 3}. (It is easy to construct similar examples when chari^ G {2, 3} as 
well.) 
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So under our last specialization, this becomes 12w(i o,i) — 12u(o,i,i). Thus the coordi- 
nates of the sole root of G in 7conv{A') — P|: are precisely [1 : — 1 : : 0]. To conclude, 
by virtue of our (compatible) choice of A' , the natural embeddings of {K*Y into Tp 
and 7conv(A') are identical. This is how we've recovered our root in Tp. 

In closing, note that in practice we would never actually compute Tiiu; s) — we 
would instead recover JF4 via rapid and sophisticated interpolation techniques, e.g., 
|Zip93| . In particular, our calculations can be sped up tremendously with suitably 



specialized code. 

2.4. The "Dense" Case. Our last example illustrates a simple fundamental case. 
Suppose E is the n-tuple ((iiA, . . . , rf„A) where A C M" is the vertex set of the 
standard n-simplex in W^ and (ij G N for all i. It is then easily verified that the n-tuple 
D := {{O , diCi} , . . . ,{0, rf„e„}) is an irreducible fill of E (cf. theorem |I|). Letting 
v4 = A, we see that our polynomial ?^ is a variant (over a general algebraically closed 
field) of the original GCP applied to an n x n system of polynomials with degrees 
(ii, . . . , (i„ | (Jan90 ]. In particular, our F — sF* has 2n s-monomials, compared to the 
n s-monomials in Canny's (/i — sx^ ^ ... , fn — sx^"). Note also that Conv(A) and 
P are homothetic and Tp = P^. Neglecting the extra s-monomials, setting di = 1 
for all i, and suitably specializing the coefficients of g, we can then recover the usual 
characteristic polynomial of a matrix. 

3. Filling 

Here we briefly recount filling and some related concepts. Some of the material 
below is covered at greater length in |[Roj94|| . The paper ||Stu94|| is also a useful 



reference but deals more with the sparse resultant than with filling. The results below 
form the basis for our combinatorial approach to perturbing degenerate polynomial 
systems. 

Let 5"^^ C M" denote the unit (n — l)-sphere centered at the origin. For any 
compact i? C M" and any w G MJ^, define B^ to be the set oi x^B where the inner- 
product x-w is minimized. (Thus B"^ is the intersection of B with its supporting 
hyperplane in the direction w.) We then define E'^ := {E^, . . . , E^) and DdE"^ := 

(DinEr,---,^nnE-). _ 

Recall that the dimension of any 5CIR", dimS, is the dimension of the smallest 
subspace of M" containing a translate of B. The following definition is fundamental 
to our development. 

Definition 4. Suppose C := (Ci,... , C„) is an n-tuple of polytopes in M" or an 
n-tuple of finite subsets o/M". We will allow any Ci to be empty and say that a 
nonempty subset J C [n] is essential for C (or C has essential subset J) <^==^ 
(0) Supp(C) D J, (1) dim(^^.g^Q) = \J\ - I, and (2) dim(^^.gj, Q) > \J'\ for all 
nonempty proper J' ^ J . 
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Equivalently, J is essential for C <^=^ the | J|-dimensional mixed volume of {Cj \ j G 
J) is and no smaller subset of J has this property. Figure |l| below shows some simple 
examples of essential subsets for C, for various C in the case n = 2. 



Ci C2 



{1},{2} 






c 



{1} 



{1,2} 



None 



Figure 1 . The essential subsets for 4 different pairs of plane polygons. 
(The segments in the third pair are meant to be parallel.) 

A basic fact about mixed volumes is that A4{E) = <^=^ E has an essential subset, 
whenever Supp(£') = [n]. However, there is an even deeper connection between filling 
and essentiality: 



Theorem 1. | Roj94 , sec. 2.5] Suppose D and E are n-tuples of finite subsets ojTP' 
such that M{E) > 0. Then D fills E ^^ for all w G >S"-S Supp(D n E"") contains 
a subset essential for E^ . D 

Remark 4. One certainly need not check infinitely many w. In fact, we need only 
check one w (just pick any inner normal) for each face of the polytope Y17=i Conv(i?j) . 

Filling is closely related to root counting for sparse polynomial systems, and this 
aspect is explored much further in Roj94 , RW96, Roj97b|| . We also point out that the 
computational complexity of finding an irreducible fill is an open question. However, 
for n<3, finding irreducible fills is quite simple and no harder (asymptotically) than 
finding a convex hull. (Using theorem |l|, this follows as a simple geometric exercise.) 
In any event, the connection between fills and polynomial system solving (not to 
mention specialized resultants) appears to be new and, we hope, provides added 
incentive to investigate filling. Also, even if finding an irreducible fill is too hard, 
this step of our toric GCP construction need only be done once for a given family 
of problems, provided E remains fixed. The situation where the monomial term 
structure of a polynomial system remains fixed once and for all, and the coefficients 
may vary many thousands of times, actually occurs frequently in many practical 
contexts such as robot control or computational geometry. 

4. ToRic Geometry and the Proof of Our Main Theorem 
Our notation is a slight variation of that used in 



Ful93|| , and is described at 
We will assume the reader to be familiar with normal 
fans of polytopes and the construction of a toric variety from a fan or a finite point 
set | Ful93|, |GKZ94 |. However, we will at least list our cast of main characters: 



greater length in |[Roj97b| 
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Definition 5. ||Roj97b|| Given any w&M.^, we will use the following notation: 

T = The algebraic torus (K*)"' 

P^ = The face of P with inner normal w 

aw = The closure of the cone generated by the inner normals of P^ 

Uw = The affine chart of Tp corresponding to the cone a^ of Fan(P) 

Lw = The dim^P"^)- dimensional subspace o/M" parallel to P"^ 

Xw = The point in U^ corresponding to the semigroup homomorphism a^ fl Z" — > 

{0, 1} mapping p \-^ S^-pfi, where 6ij denotes the Kronecker delta 

Oyj = The T -orbit of x^ = The T -orbit corresponding to RellntP'" 

^p{Q) = The T -invariant Weil divisor of Tp corresponding to a polytope Q. 

Div(/) = The Weil divisor of Tp defined by a rational function f on (K*)^ 

Vp{f,Q) = Div(/) -\- £p{Q) = The toric effective divisor ofTp corresponding to {f,Q) 

Vpi^F, V) = The (nonnegative) cycle in the Chow ring ofTp defined by nj=i 1^p{fii Pi), when- 
ever V = {Pi, .. . ,Pk) 

We will say that a polytope P is compatible with Q iff every cone of Fan{Q) is 
a union of cones of Fan(P) ||Kho77 , i'ulQSI , |Roj97b |. In particular, whenever F is an 



k X n polynomial system with support contained in E, we will define the zero set 
of F in Tp to be the toric cycle Vp{F,V), where P:=(Conv(Ei), . . . , Conv(Efc)). 

The following result will provide some necessary geometric intuition for specializing 
resultants. 



Vanishing Theorem for Resultants. [ [Roi97d|] Suppose fi is a polynomial over K 



with support contained in Ei C Z" for all z G [n + 1] . Then, provided 

M{Ei,... , Ei_i, Ei^i, . . . ,En+i)>0 
for some iG [n + 1], 

71+1 

Res^(/i,... ,/„+i) = 0^^ f|Pp(/„Conv(i?,))^0, 

i=l 
^ n+1 

where E:={Ei, . . . ,i?„+i) and P'-=Y1 Conv(_Ej). D 

j=i 

Remark 5. This result provides a geometric analogue, over a general algebraically 
closed field, of the product formula for the sparse resultant [pS93|| . 

We will also make frequent use of the natural correspondence between the face 
interiors {RellntP"'} and the T-orbits {O^} | |KSZ92| , fW{ , ^KZg. The following 



lemma gives a more explicit algebraic analogy between the faces of P and the affine 
charts of Tp. 



Lemma 1. |P^oj97b| , Sections 4.2-5.1] Suppose F is a k x n polynomial system over 



K with support contained in a k-tuple of integral polytopes V := {Pi,... ,Pk) in 
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M". Assume further that P is a rational polytope in R". Then the defining ideal 
m K[x^ \e e alnZ""] of U^ n Vp{F,V) is {x^'fi \ for all i G [A;] and b, e Z" such 
that hi + Pi G al) . U 

Lifting (or projecting) from one toric variety to another is an important fundamen- 
tal ideal we will also use. The following lemma follows directly from the development 
of iF\Il93l . 



Lemma 2. Suppose P C M" is an n-dimensional rational polytope, and A is either 
a nonempty finite subset of 7/'' or a rational polytope in W^ . Assume further that P 
is compatible with Conv(74). Then there is a natural (surjective) proper morphism 
(f : Tp ->^ Ta- In particular, following the notation of this section, (f{Vp{F,V)) = 
Va{F,V), where the latter cycle is the image o/ T'conv(yi)(-^, "P) under the natural 
proper morphism from Tconv{A) to Ta- 

Remark 6. Recall that Ta can be defined as the image of of Tp under the map Lp 
from our Main Theorem. So, with this understanding, there is no ambiguity between 
our first and second if. 

To conclude our background, we will need the following lemma implying that F* 
is sufficiently generic in a useful sense. 

Lemma 3. Suppose D is an irreducible fill of some n-tuple E. Then for any point v 
lying in any Di, there exists a w gW^\{0} such that {i} is the unique essential subset 
of D^ and Df = {v}. In particular, following the notation of our Main Theorem, F* 
has exactly Ai{V) roots (counting multiplicities) in {K*)^ . D 

This lemma follows easily from the techniques of [[Roj94|| , particularly section 2.5. 

4.1. The Proof of Our Main Theorem. We first note that the well known results 

on the degree of Res5(/i, . . . , /n+i) with respect to the coefficients of different fi 

Stu94|| remain true over any algebraically closed field. This follows easily from the 



formulation of the resultant for a collection of invertible sheafs on a projective variety 
|pKZ94|| . In particular, the degree of 7i as a polynomial in s should be 

n 

2_^M.{Ei, ■ ■ ■ ,Ei^i, -Ej+i, . . . ,En,A). 

i=\ 

Also each coefficient of ?i(s) should be a homogeneous polynomial (in the Wg) of 
degree M.{E'). These two assertions of course include the opening statement of our 
Main Theorem (on the degree and homogeneity of Ta)^ but they will follow only 
upon showing that 7i is not identically 0. 

To see this, note that lemma |l| and the Vanishing Theorem for Resultants readily 
imply that the coefficient of the highest power of s in 7i is precisely Res(^^yi)(F*, g). 
(Simply check the zero set oi F — sF* in Tp at s = oo.) By lemma ^ and the Vanishing 
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Theorem once more, we see that this polynomial in the Ue is not identically 0. So 
Ti^O and we've finished the simplest part of our proof. 

Part (0) of our Main Theorem follows similarly: One need only consider the un- 
specialized resultant polynomial Iies(^E,A){F,g) and observe the terms of degree 
in s as we specialize coefficients to obtain F — sF*. The statement on the vanishing 
of Ch^ then follows easily from lemma ^ (since P is compatible with Conv{A)) and 
the Vanishing Theorem: We obtain that (p{Z) is positive-dimensional iff Ch^ has 
infinitely many distinct divisors of the form YleeA'-^eiJ-e- In particular, corollary |^ 
follows as a special case of (0) since P compatible with Conv{A) =^ Tp = Tp. Note 
that (1), (2), and (3) also follow almost trivially, provided Ch^ is not identically 0. 

To properly handle the cases of (1), (2), and (3) where we are actually working with 
a non-trivial toric perturbation, let us first construct two important toric cycles: Let 
Z be the zero set of F - sF* in Tp x P^^ and Zl the zero set of n{u; s) in pl^'"^ x P}^.. 
Then it is easily observed that Z = Z n {Tp x {0}). Also, it can be shown |[Roj97b| , 



Section 5.1] that Z contains an algebraic curve C (possibly reducible), with surjective 
projection onto the second factor oi Tp x F]^, obeying the following property: C fl 
[Tp X {0}) contains the zero-dimensional part of Z, and consists of exactly A4{E) 
points (counting multiplicities). (In fact, Zn (Tp x {sq}) = Cn {Tp x {sq}) for almost 
all So € P/f.) Furthermore, by slightly modifying step (i) of the proof of the toric 
variety version of Bernshtein's theorem ||Roj97b| , Section 5.1], one can show that C 



intersects every positive-dimensional irreducible component of Z. (One also needs to 
use the definition of intersection multiplicity of an irreducible component W as the 
number of curve branches intersecting VF in a 1-parameter deformation.) 

The proof of the rest of our main theorem will reduce to establishing a precise 
correspondence between the factors of J-'a and the points of (fi{C fl (Tp x {0})). To 
complete this connection, we need only observe that Z^ is a very special kind of 
hypersurface, closely related to Z. 

In particular, if k is the least power of s in H, observe that Z^ and the zero set of ^ 

in Pj^ X F]^ differ only by the presence of the hyperplane P^^ ^ {0}- The second 
zero set does not contain this hyperplane, so let's call the second zero set Z^. Then by 
lemmata |l] and H, and the Vanishing Theorem for Resultants, dim[Z' fl [Tp x {sq})] = 
implies the following equivalence: H(H^(^i^y, So) = <^=^ (eZ H (TpX{so}), where Hp 
is the hyperplane dual to the point p.f\ Now dim[Zn (Tp x {oo})] = by construction. 
Hence dim[Zn (7p x {sq})] =0 for almost all SqEF]^, by Main Theorem 2 of |[Roj97b| . 



So Tp^ (C) is an open subset of Z^ , where we define 4^^ (C) :={(y, Sq) \ yGH^i^Q ; ^G 
C n (Tp X {so}) ; So gP]^}. Therefore, since tp is a proper map, ^ must vanish on 



^So iip:=[pe I eeA]epl^l ' then Hp:={[ye \ eeA]eP^I ' \ EeeAPeVe^^}- 
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all of '0^(C). In particular, 

Cecn(r^x{o}) VeeA 

where a&K* , [c^_e | eeA] ■.= ip[C,), and the product counts intersection multiplicities. 

Continuing our main proof, (1), (2), and (3) follow immediately from our last 
formula and our preceding observations. D 

Note that our algebraic proof avoids the use of limiting arguments that were present 
in | ](Jan9(J| | . Thus our result is equally valid when K has positive characteristic. 
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